Dynamics and planet formation in/around 
binaries 



F. Marzari^ P. Thebault^'^ S. Kortenkamp'' and H. SchoU'^ 



4 



Dipartimento di Fisica, Via Marzolo 8, 35131 Padova, Italy marzari@pd.infn.it 
Observatoire de Paris, Section de Meudon, F-92195 Meudon Principal Cedex, 
France philippe .thebault@obspm.fr 

Stockholm Observatory, Albanova Universitetcentrum, SE-10691 Stockholm, 
Sweden thebaultSastro .su.se 

Planetary Science Institute, Tucson, Arizona kortenka(§psi . edu 
* Observatoire de la Cote d'Azur, Nice, France scholl@obs-nice.fr 

Summary. 

We study to which extent planetesimal accretion is affected by the perturb- 
ing presence of a compagnon star. We concentrate on one crucial parameter: the 
distribution of encounter velocities within the planetesimal swarm. We numerically 
explore the evolution of this parameter taking into account the secular perturbations 
of the binary and friction due to the very likely presence of gas in the disk. < Av > 
maps are derived, for planetesimals of different sizes, for a total of 120 binary con- 
figurations (eccentricity eb and separation ab) . We identify for each case 3 different 
accreting behaviours. 1) In regions where no significant dV increase is observed, 
"classical" (i.e., single-star) accretion is possible. 2) In regions where dV > Ve.ro, the 
threshold velocity above which all impacts are eroding, no accretion is possible and 
planet growth is stopped. 3) In between these 2 limiting behaviours, a large fraction 
of binary configurations leads to significant dV increase, but still below the erosion 
threshold. In this intermediate case, planetesimal growth can occur, but proceeds 
slower than in the single-star case, possibly following the so-called type II runaway 
groth mode. 



1 Introduction 

1.1 the "standard" planet formation scenario 

A fundamental step in the "standard" scenario for the formation of planets, 
either terrestrial or gaseous giant, is the coUisional accumulation of kilometer- 
sized planetesimals in successively larger bodies. It is widely admitted that 
planetesimals emerge from the coagulation of silicate and ice grains in the 
midplane of protoplanetary disks on timescales of the order of 10^ yrs, even 
if the details of the formation process are not well understood at present. 
Colliding with each other, planetesimals grow into large planetary embryos 
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(Lunar- to Mars-size) through a phase of "runaway" and "oHgarchic" growth 
on a timcscalc of the order of 10"* — 10"'' years (Grecnberg et al., 1978; Wetherill 
& Stewart, 1989; Barge & Pellat, 1993; Lissauer, 1993; Kokubo & Ida, 1998, 
2000; Rafikov, 2003, 2004). This is a robust finding corroborated by distinct 
numerical methods (see review by Kortenkamp et al., 2000). For terrestrial 
planets, the formation process reaches completion after a final phase of giant 
impacts and mutual accretion of planetary embryos resulting in 10^'^ — 10^^ 
gram full size bodies in about 10^ — 10* years on stable separated orbits 
(Chambers and Wetherill, 1998). In the case of giant planets, protoplanets 
accumulate into a solid core that, once sufficiently massive, begins to ac- 
crete a gaseous envelope. Initially, this envelope is in hydrostatic equilibrium 
sustained by the luminosity provided by the accreting planetesimals and is 
increasing at a slow rate compared to that of the solid core. However, at a 
critical core mass (in between 5 to 15 M(q) a phase of runaway gas accretion 
occurs, the planet mass increases rapidly and it reaches its final value on a 
timescale estimated to be only a few 10^ yr. For a given density of solids and 
gas in the disk, hydrodynamical calculations can determine for which mass 
of the core the envelope is able to remain static or collapse, and the rate at 
which the mass flows onto the planet (see e.g. Wuchterl et al., 2000). At the 
end of the gas rapid infall, the planet evolves into the isolated stage during 
which it cannot grow any further because of the exaustion of the gas sup- 
ply either because of its removal by effects of the star or as a result of gap 
formation around the planet orbit triggered by disk-planet tidal interaction 
(Papaloizou and Lin, 1984). The planet finally cools down and contracts to 
its present size. This scenario for giant planet formation, called core accretion 
(Pollack et al, 1996), is based on the interior models for Jupiter and Saturn 
which strongly suggest the presence of a dense central core surrounded by an 
envelope composed mainly by hydrogen and helium for both planets. There 
are large uncertainties on the mass of the core and detailed calculations per- 
formed by (e.g. Guillot, 1999) predict un upper limit of 10 Mq for Jupiter's 
core while for Saturn the core mass ranges from 6 to 17 M^. 

It is believed that extrasolar planets, detected in orbit around an increasing 
number of nearby stars may have formed either by core accretion or by disk 
instability (e.g. Boss, 1997). While core accretion requires a few million years 
to form a gas giant planet, in the competing model self-gravitating clumps 
formed by disk instability will contract to planetary densities in times of a 
few hundred years. The two formation mechanisms are both viable and, at 
the moment, there is no reason to definitively exclude one in favor of the 
other. In this paper, however, we will concentrate on planetesimal accretion 
as a common mechanism for forming either terrestrial planets or the core of 
gaseous giants. 

In a given protoplanetary disk, the potentiality of a planetesimal swarm 
to accumulate into larger bodies is measured by two critical parameters. 
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• The surface density of solid material. The concentration of condensible 

material determines the final mass in the case of a terrestrial planet, while 
it must lead to a reasonable formation time for the core of a giant planet 
that must be compatible with the lifetime of the gaseous disk. Disk ages 
are estimated to range from 0.1 to 10 Myr (Strom et al.. 1993; Haiscli 
et al., 2001; Chen and Kamp, 2004), while 3 Myr is the age at which 
half of the stars show evidence of disks. Simulations of the accumulation 
of jovian planets (Bodenheimer and Lin, 2002) show that if the surface 
density of solids in the disk is assumed to be at least 3-times larger than 
that of the minimum-mass solar nebula (Weidenschilling, 1977) then the 
formation time for Jupiter mass planets is within 2 3 Myr at 5 AU from the 
star. If planetary migration is included in the simulations, the timescale 
for forming Jupiter is reduced to 1-2 Myr even with solid surface density 
closer to the minimum mass solar nebula (Rice and Armitage, 2003). 

• The relative encounter velocities between planetesimals. 

The random encounter velocities in a planetesimal swarm determine 
whether mutual impacts will result in accretion or on the contrary in 
cratering or fragmentation. The key mechanisms controlling the impact 
velocities are: 1) mutual gravitational encounters between planetesimals, 
2) physical collisions, 3) gas drag, in presence of the gaseous component 
of the disk. Several numerical studies have shown that in an unperturbed 
swarm of kilometer-sized planetesimals around a single star the average 
random encounter velocity (Av) is low enough to favor the coUisional ac- 
cumulation of bodies (e.g. Greenberg et al., 1978). At the early stages of 
accretion a low (Av) is also a fundamental requirement for permitting a 
phase of rapid " runaway" growth which significantly shortens the timescale 
for the formation of planetary embryos. Indeed, under this condition the 
ratio Vesc{R) I {Av) between the escape speed of the largest bodies and the 
average impact velocities is >> 1, thus making their gravitational cross- 
section quite larger than the geometrical cross-section. This gravitational 
focusing enhances the growth rate of the larger planetesimals which grow 
much faster than the rest of the population (e.g. Lissauer, 1993). A direct 
consequence of this crucial {Av) dependence is that planetesimal growth is 
very sensitive to any increase of the encounter velocities. It could thus be 
significantly slowed down or even stopped if the gravitational influence of 
a massive body external to the swarm, like a a pre-existing giant planets 
(e.g., Jupiter in our solar system) or a companion star in a binary sys- 
tem is able to excite large impact velocities. Thcbault et al. (2002) have 
shown that the formation of terrestrial planets in those exoplanetary sys- 
tems hosting a giant planet on an external orbit may have been indeed 
halted. They find that in both the e Eridani and 47 Uma planetary sys- 
tems, the relative velocity between planetesimals in the terrestrial planet 
region may have quickly exceed the critical threshold value v^ro beyond 
which erosion dominates before the completion of the runaway growth if 
the massive outer planet rapidly reached its present mass. Of course, the 
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threshold velocity that preferentially result in disruption rather than ac- 
cretion is a sensitive function of the internal strength of planetesimals, 
their bulk density and the fraction of impact energy partitioned into the 
fragments after the impact. These physical parameters strongly depend on 
the initial radial density and temperature profile in the disk midplane. 

1.2 planets in binaries, planet formation in binaries 

In this paper we will focus on how the presence of a companion star affects the 
planetesimal accretion process either in the formation of terrestrial planets or 
the core of a giant planet. While the existence of extrasolar Earth like planets 
has yet to be assessed, some of the giant exoplanets (up to about 40 so far) have 
been detected in the so-called S-type orbits that encircle one component of a 
binary system (e.g. Raghavan et al., 2006; Desidera & Barbieri, 2007). Most 
of these are widely separated pairs where the stellar companion probably does 
not have a significant influence on the planetesimal accretion around the other 
star. However, in three of the systems, •jCephei GL 86, and HD41004, the 
companion is within 20-25 AU from the star hosting the planet. The ^Cephei 
system is made of a central star with an estimated mass of 1.40 ± O.I2M0, a 
planet orbiting at 2.13±0.05 AU with a minimum mass of 1.6±0.13 Mj and a 
companion star on an eccentric orbit (e = 0.411 ±0.0063) with semimajor axis 
of 20.23±0.64 AU and mass of about O.4M0 (Hatzes et al, 2003; Torres, 2007; 
Neuhauser et al., 2007). GL 86 is a dwarf star somewhat less massive than the 
Sun (about O.7M0, e.g. Santos et al., 2004) with a white dwarf companion 
of O.55M0 (Mugrauer and Neuhauser, 2005) on an orbit having a semimajor 
axis of 18.42 AU and an eccentricity of e=0.3974 (Lagrange et al., 2006). A 
massive jovian planet {Msin{i) = 3.91Mj) moves around the primary star on 
an almost circular orbit (e = 0.0416) with a semimajor axis of a = 0.113 AU 
(Eggenberger et al., 2003). HD41004 A is a O.7M0 star with a dwarf M4V 
companion orbiting at ~ 22^417. The primary star has a Msin{i) = 2.3Mj 
planetary companion on a very eccentric e=0.39 orbit at a = 1.31 AU (Zucker 
et al., 2004; Raghavan et al., 2006). Note that this system is a hierarchical 
quadruple, since the stellar companion is itself orbited at a = 0.016 AU by a 
brown dwarf of minimum mass 18.4Mj (Zucker et al., 2004). 

The observed characteristics of these systems, e.g., the proximity of the 
companion star and its eccentric orbit, strongly suggest that the presence of 
the secondary star must have had an influence on the formation of the detected 
planets. The gravitational pull of the companion star acts on the (Av) of the 
planetesimal disk around the primary star. As previously mentioned, this can 
modify the course of the accretion process which strongly depends on the ratio 
VescI {'Av) , where v^ac is the surface escape velocity of the biggest accreting 
bodies. Timing is an important issue since the perturbations will be effective in 
altering the planetesimal formation process only if tgec, the timescale required 
by the the secular perturbations of the star to induce a significant Av increase, 
is shorter than tgrow, the timescale for the runaway/oligarchic formation of 
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embryos which is typically of the order of 10* to lO^yrs. We can outline three 
possible different evolutions for a perturbed planetesimal swarm: 

• 1) If (Av) remains smaller than Vesc during the growth phase then Plane- 
tary formation proceeds almost unaffected by the companion perturbations 
. The "classical" planetesimal accumulation scenario holds. 

• 2) For values of (Av) larger than Vesc but still smaller than Veroj the 
threshold velocity for erosion to dominate over accretion, planetary accre- 
tion will still be possible but runaway growth will either not occur or will 
occur only after large enough bodies have formed such that Vggc > {Av) 
(the so-called Type II runaway growth described in Section 6 below). 

• 3) If (Av) is increased beyond Vero, in the majority of collisions cratering 
and fragmentation will overcome accretion and the planetesimal popula- 
tion will be slowly ground down to dust, failing to form a full size planet. 

In the following, wc will focus on the determination of (Av) for differ- 
ent parameters of the binary star system with the goal of evaluating when 
planetesimal accretion is possible. We have performed a series of determin- 
istic numerical simulations, following the dynamical evolution of a swarm of 
test planetesimals, under the influence of the companion star's gravitational 
pull and friction by the gas in which the planetesimals are imbedded. Crucial 
parameters, such as the companion star's mass mi,, semi-major axis at, and 
orbital eccentricity Cb are explored as free parameters. We use a code initially 
developed for the study of planetesimal systems perturbed by giant planets 
(Thcbault & Brahic, 1998; Thebault & Beust, 2001; Thcbault ct al., 2002) 
and adapted to the circumprimary case (Marzari & SchoU, 2000; Thebault et 
al., 2004, 2006). This 3-D code has a close encounter tracking algorithm which 
enables precise determination of the Av within the system. 

In a first step we present a detailed study of the planetesimal secular dy- 
namics in presence of a companion star (section 2). We will then proceed 
by introducing the effects of gas drag and by estimating how the combined 
action of the friction force with the secular perturbations affect the value of 
{Av) (section 3). We then compare the derived {Av) with an an approximate 
estimate of Uero to derive a yes-or-no criterion for planetary formation as a 
function of the binary orbital parameters (section 4). Finally, in section 5, 
wc present and describe one possible accretion-growth mode for planetesi- 
mals in perturbed binary star systems, the so-called Type II runaway growth 
mechanism. 



2 Planetesimal dynamics in a binary star system: the 
secular approximation 

Hereinafter, we will assume that the orbit of the secondary star is coplanar 
to the planetesimal disk so that we can neglect the inclination of the bodies. 



6 



F. Marzari P. Thebault S. Kortenkamp and H. Scholl 



This is a reasonable assumption when the star is close to the primary because 
of the short relaxation time of the disk into the plane of the binary orbit. 

Heppenheimer (1978) developed a simplified theory for the evolution of the 
planetesimal eccentricity with time based on an expansion of the disturbing 
function in a power series of the ratio of the semi-major axes a/af, where a 
is the semimajor axis of the planetesimal and Ub that of the companion star. 
The approximation holds if a is small compared to Ub and it stays within 
the critical semimajor axis Oc for dynamical stability derived using direct 
numerical integrations by Wiegert and Holman (1997). By truncating the 
disturbing function of the secondary star at the second order in the eccentricity 
of the planetesimal, the seciilar equations simplify and lead to the following 
expressions for the forced eccentricity ep and the frequency u of circulation 
of the perihelion w. 

5 a 66 
4 afe 1 - eg 

3 1 

4(1- e2)3/2 al 

with 

tan rot = ^ ' , 3 

^ ^ " 1 - cos(2ut) ^ ' 

However, Thebault et al. (2006) compared the predictions of the above 
equations with direct numerical simulations and found significant discrepan- 
cies, especially for the oscillation frequency u. They empirically derived a 
revised expression for the frequency u which reads: 



3 1 a3/2 



1 + 32. 



1 - ef \ab 



(4) 



The discrepancy between the numerical and analytical estimate of u is reduced 
from about 70%, in most of the cases with high eb and e, to less than 5%. 

To illustrate the dynamical evolution of a planetesimal swarm under the 
secular perturbations of a companion star we have performed a "pedagogical" 
numerical simulation where we integrate the orbits of massless test bodies, 
initially on circular oribts, for about 10^ yrs. In Fig.l we plot the location 
of each particle in the (a,e) plane at different times. The eccentricity of each 
single body oscillates with a period u that depends on the semimajor axis 
of the body. If we freeze the dynamical system at subsequent times, the slow 
dephasing of the eccentricity of nearby orbits will translate into a wavy pattern 
whose density grows with time. 
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Fig. 1. Evolution of a test particle populatiou perturbed by a stellar companion hav- 
ing approximately the characteristics of 7 Cephei: with mt = 0.25, at = 20 AU and 
eb = 0.3. The dotted line represents the distribution of average encounter velocities 
between the bodies of the system. 



3 Effects of the secular perturbations on (Av) 

If the orientations of the perihelia were fully randomized in the planetesimal 
swarm, the effect of the secular perturbations of the companion star on {Av) 
would be easy to derive through the standard relation (Lissauer & Stewart, 
1993): 

(Av) ~ {5/A)'/^eFVK (5) 

where vk is the local Keplerian velocity. According to this equation, even 
very low values of forced eccentricity would halt accretion. In the 7 Cephei 
binary system, for example, {Av) would be approximately 1.8 x 10"' m.s^-'^ at 
2 AU, the distance from the star where the planet has been found. This large 
relative velocity would certainly prevent mutual accretion in any swarm of 
kilometer-sized planetcsimals (which have escape velocities of the order of a 
few m.s~^). In addition, a proper eccentricity component Cp has to be added 
to the planetesimal ecentricities, causing oscillations of the eccentricity around 
ep with an amplitude equal to Cp, which leads to even higher impact velocities. 
If this was the scenario in the early phases of planetesimal accumulation in 
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close binary systems, secular perturbations pose a serious threat to planetary 
formation. However, the predictions of Equ.5 strongly depend on the assump- 
tion that planetesimal trajectories have a randomized distribution of orbital 
angles. This is not true at the beginning of the dynamical evolution, because 
secular perturbations force of a strong phasing of the orbits. The large eccen- 
tricity oscillations induced by the companion star are in phase and associated 
to a very efficient periastron alignment between neighboring orbits. This con- 
figuration leads to planetesimal encounters on almost tangential trajectories 
and, notwithstanding the very high values of forced eccentricity, to low and 
accretion friendly relative velocities, dominated by the Keplerian shear. 

At subsequent times, however, the dependence of the precession frequency 
u of CC7 on a^/^Aa causes a progressive dephasing of nearby orbits. At some 
point neighbouring orbits eventually cross (see the detailed discussion in 
Thebault ct al. (2006)), leading to a sudden and very sharp increase of 
the encounter velocities. This is illustrated in Fig.l where we have super- 
imposed to the (e,a) graph the corresponding evolution of the average rela- 
tive velocity (Av). The transition between the inner regions, where the or- 
bital phasing is still strong preventing orbital crossing, and the outer regions 
where orbits already cross is very abrupt: (Av) changes from a few m.s~^ 
to about 1km. s""'^, an accretion-inhibiting value for kilometer sized planetes- 
imals, within ~ 0.1 AU. Thebault et al. (2006) have derived an analytical 
estimate of the location in semimajor axis across where the orbital crossing 
occurs as a function of time: 



Cia, 



1/2 



+ C2 

ab J \ ab 



, 4 
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where a^, e{, and are the binary semimajor axis, eccentricity and mass 
ratio, respectively, and: 

^ 457r et 1 „ 224 mt 

Conversely, for a given location a in the planetesimal swarm, the time tcr 
required by the wavefront of high impact speeds to reach across is: 



« no ^ in2(l -e^)^ / "^b \'^-^ /across \-^-^ 

W = 6.98xl0^-^(^— j (^) (^) yr. (8) 

These analytical formulae hold for at > 10 AU and for 0.05 < eb < 0.8. In 
Fig. 2 wc; show the time required by the high speed regime to reach 1 AU for 
for different values of Of, and ef, and for nib = 0.5. Note that for low Cf, and 
large ab, tcr is longer than 10^ yr, giving a chance to runaway growth to build 
up large protoplanets before an erosion regime is established. 
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Fig. 2. Value of the minimum companion semi-major a^xis, ab.cri, leading to orbital 
crossing of planetesimals at 1 AU, for different values of the crossing time tcr, as a 
function of the companion eccentricity. The companion's mass is fixed, with m.b = 
0.5. The dotted lines represent constant values of (Av) at orbital crossing, as given 
by Equ.5 

An implicit assumption of this scenario is that the initial encounter ve- 
locities in the planetesimal swarm are very low, i.e. that the initial orbits 
are either circular or fully phased. This assumption is commonly adopted 
when addressing planetesimal accretion in binaries (e.g. Heppenhcimer, 1978; 
Whitmire et al., 1998; Thebault et al., 2004), but should however be handled 
with care. Its validity depends on how and how fast planetesimals decouple 
from the gas of the protoplanetary disk to follow their own Keplerian orbits 
around the main star. If the decoupling occurs almost at the same time for all 
the planetesimals and if they are at that moment on almost circular orbits, 
then their proper eccentricities Cp's are approximately all equal to ep and 
the w's are grouped around with consequent low relative velocities. Basi- 
cally, it all goes down to how the "initial" kilometer-sized planetesimals are 
being formed. Unfortunately, so far there are no clues on how planetesimals 
form in a disk perturbed by the gravity of the secondary star and how these 
perturbations may affect their trajectories as they emerge from the gas of the 
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disk. As a matter of fact, there is not even full agreement on how planetesimals 
form around "normal" quiet single stars. This very difficult exceeds by far the 
scope of the present, but there are schematically two competing models for 
planetesimal formation in a quiescent protoplanetary disk: gravitational insta- 
bility (e.g. Goldrcich & Ward, 1973; Youdin & Shu, 2002; Youdin & Chiang, 
2004) , where kilometer-sized bodies directly form from small solids in a dense 
unstable solid grain layer in the midplane of the protoplanetary disk, and 
collisional-chcmical sticking (e.g. Wcidcnschilling, 1980; Dominik & Tielens, 
1997; DuUemond & Dominik, 2005), in which planetesimals are the outcome 
of a progressive mutual grain sticking. If gravitational instability works out 
also for disks in binary systems, then the planetesimal formation timescalc is 
likely to be negligible compared to the typical timescales for runaway growth 
and for the onset of fully dephased secular perturbations by the companion 
star. In this case, it is reasonable to expect an initial orbital alignment within 
the swarm. The decoupling from the gas occurs at a slower rate in the chemical 
sticking scenario where planetesimals grow progressively. They might begin 
to feel the companion star perturbations when they are not yet fully detached 
from gas streamlines. However, numerical simulations indicate that the whole 
growth process from grains to kilometer-sized bodies might not exceed a few 
10^ years ( Wcidcnschilling, 2000), i.e. still shorter than both runaway growth 
and secular perturbations timescales. 



4 Role of gas drag 

So far, all effects other than the gravitational perturbations of the compan- 
ion star have been neglected. However, in the standard planetary formation 
scenario, the initial stages of planetesimal accumulation are believed to take 
place in presence of significant amounts of primordial gas. Marzari and Scholl 
(2000) have shown that gas drag damps planetesimal eccentricities and also 
restores a strong perihelion alignment for equal size bodies. If the density of 
the gas is large enough all the tus of planetesimals halt their circulation and 
tends towards a fixed value Wai equal to 270°. This tends to cancel out the 
large eccentricity oscillations forced by the companion. In Fig. 3 we show this 
phenomenon for 5 km planetesimals embedded in a protostellar disk similar to 
the minimum mass solar nebula and perturbed by an outer star. The parame- 
ters of the binary system recall those of a Centauri. In this simulation and in 
all the following ones, the drag force is modelled in laminar gas approximation 
as 

r = -KVrel-Vrel, (9) 

where K is the drag parameter given hyK = 3 10 (Kary et al., 1993) where 
s is the radius of the planetesimal, ppl its mass density, pg the gas density of 
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the protoplanetary disk and a dimensionless drag coefficient related to the 
shape of the body 0.4 for spherical bodies) ^ 

The orbits of the planetesimals are all aligned for an indefinite time up to 
about 2 AU from the star and the {Av) is of the order of a few meters per 
second. Under this condition, planetesimal collisions lead to accretion all the 
time. 

Unfortunately, this alignment depends on the balancing between the drag 
force and the strength of the secular perturbations. For a lower gas density 
or for a larger size of the planetesimals the perihelia moves away from vjai 
towards larger values and the alignment is perturbed by a wavy pattern that 
forewarns the onset of a regime dominated by secular perturbations. This 
behaviour is illustrated in Fig. 4 for a bimodal planetesimal population char- 
acterized by two different sizes. 

The different alignment of planetesimals according to their size has a cru- 
cial impact on the accretion process. While equal size planetesimals keep a low 
{Av) since their orbits are oriented towards the same angle, when different 
size planetesimals encounter each other the relative velocity is large most of 
the time because of the forced unpaired orbital alignment. The high eccen- 
tricities induced by the secular perturbations introduce a large non-tangential 
component in the relative velocity that easily overcomes the Keplerian shear 
and cause high velocity impacts. In Tables I and II we report the average 
encounter velocities < Z\w(i?i, R2) > between planetesimals of size Ri and R2 
for two different binary star orbital configurations. The assumed density of 
the gas drag corresponds to that of the Minimum Mass Solar Nebulae derived 
by Hayashi (1981). The velocity values have been obtained with numerical 
integrations of the planetesimal trajectories over 2 x lO'^years. For the case 

^ Our gas drag model is a simplified one where the gas disk is assumed to be fully 
axisymmetric and follows a classical Hayashi (1981) power law distribution. It 
is however more than likely that in reality the gas disk should depart from this 
simplified view because it would also "fool" the companion star's perturbations. 
Several numerical studies have investigated the complex behaviour of gaseous 
disks in binary systems. They all show that pronounced spiral structures rapidly 
form within the disk (e.g. Artymowicz and Lubow, 1994; Savonije et al., 1994) 
and that gas streamlines exhibit radial velocities. To follow the dynamical be- 
haviour of planetesimals in such non-axisymmetric gas profiles would require a 
study of the coupled evolution of both gas and planetesimal populations, have 
to rely on hydro-code modeling of the gas in addition to N-body type models 
for the planetesimals. Such an all-encompassing gas-|-planetesimals modeling is 
clearly the next step in binary disk studios and have already been started by 
several teams. It is interesting to note that preliminary results seem to show that 
planetesimal behaviours in systems with "realisitic" gas disk modeling do not 
seem to drastically depart from the behaviour in the axisymmetric case. There is 
in particular no phase alignment between eccentric planetesimal orbits and gas 
streamlines, so that gas friction on planetesimals is still very high (Paardekooper, 
private communication) . 




Fig. 3. Pcriastron and eccentricity distributions of a swarm of 5 km diameter plan- 
etesimals surrounding the main star of the Alpha Centauri system. After 10* yrs 
from the beginning of the simulation the planetesimal orbits are strongly aligned be- 
cause of the coupling between gas drag and secular perturbations of the companion 
star. 



displayed in Table I, it is worth noticing that, in absence of the gaseous com- 
ponent, the parameters of the binary system would lead to orbital crossing of 
nearby planetesimal orbits in less than ~ 5 x lO'^years (see Fig. 2). The swarm 
would stop growing and erosion would take place. When we include the effects 
of gas drag the {Av) remains always low for equal size planctcsimals. How- 
ever, for bodies of different sizes, the different phasing caused by the coupling 
between gas drag and secular perturbations pump up the relative velocities. 
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Fig. 4. Example gas drag run. Snapshots, at i = 3 x lO^yrs, of the (e, a) and 
{w — zub,a) distributions for 2 planetesimal populations of size Ri = 1km and 
J?2 = 5 km. tu — Wb is the difference, in angular degrees, between the particles and 
the companion star's longitude of periastron. The companion star orbital parameters 
are: at = 10 AU, Cb = 0.3, mb = 0.5 



Table 1. Average encounter velocities in m.s^^, at 1 AU from the primary, within 
a population of "small" planetesimals 1 < i? < 10 km for a gas drag simulation 
with the companion star parameters nib = 0.5, at = 10 AU and Cb = 0.3 nib = 0.5, 
ttj, = 20 AU and Cb = 0.4. Avm,R2 values are averaged over the time interval < 
t < 2 X lO^yrs. Initial starting encounter velocities are such as Avo ^ lOm.s."^. 
Av values in bold correspond to accreting impacts for all tested collision outcome 
prescriptions. Underlined values are those for which we obtain different accretion 
vs. erosion balance depending on the tested prescription. Values in classical roman 
characters correspond to cases for which all tested models agree on a net erosive 
outcome. 



Sizes (km) 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 



7 8 9 10 



10 154 233 

172 10 94 
238 84 11 
289 144 63 
325 188 103 
373 228 144 
400 261 182 
428 298 212 
450 310 238 
453 338 263 



285 327 
133 187 
54 99 



12 

43 
83 
113 
147 



40 
12 

32 

68 
98 



168 123 
196 152 



360 391 
223 262 
137 177 
80 115 
32 70 
11 36 
36 12 
56 35 
83 48 
107 73 



426 452 458 
287 316 334 
200 230 254 
149 171 198 
100 122 154 
56 84 104 
48 
36 
13 
31 



35 
12 
36 

48 



76 
45 
31 
13 



14 
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Table 2. Same as Table 1, but with companion star parameters mt, = 0.5, 05 
20 AU and et, = 0.4. 



Sizes (km^ 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



8 9 10 



11 127 
126 10 
200 91 

258 146 
301 192 
339 232 
361 262 
395 295 
425 320 
446 346 



204 255 298 
84 139 185 
11 68 108 
61 9 48 
113 54 12 
152 99 43 
181 126 77 
219 159 112 
246 190 136 
266 208 163 



342 368 
227 258 
158 186 
88 120 
44 75 
10 31 
28 13 
68 26 
92 55 
122 82 



390 417 442 
290 317 340 
218 246 272 
154 186 209 
111 136 164 
66 92 119 
26 56 
11 28 
25 11 
49 22 



87 
48 
23 
12 



5 Dependence of the accretion process on the binciry 
pcLrameters. 

Once we determine the {Avi^n^^n^)) for any pair of planetesimal size we must 
compare it with the threshold velocity Vero for that pair of projectile and target 
to test whether either erosion or accretion will be the result of the collision. 
As pointed out in Thebault et al. (2006), in literature one can find many 
prescriptions of the coUisional outcome between two planetesimals given their 
impact velocity. Adopting the approach outlined in Thebault et al. (2006), 
three different values of computed according to Holsapple (1994), 

Marzari et al. (1995) and Benz & Asphaug (1999). Erosion (or accumulation) 
is assumed to occur when ^j^)) is higher (or lower) than all the three 

values of Vgro- In the intermediate cases we prefer not to draw any definitive 
conclusion letting them as undefined cases. In Tabs.l&2 a different text font 
is used to indicate when the computed (Avf^n^^n^^) leads to accretion, erosion, 
or to an undefined situation. By inspecting the tables, it is interesting to note 
that anytime we depart from the diagonal stripe Ri = R2 there is a sharp 
increase of (Av), almost always exceeding the threshold value for accretion, 
at least for the binary parameters adopted in the tables. 

The knowledge of the relative velocities between different size planetesi- 
mals is the first step towards fixing constraints on the dynamical environment 
allowing planetary formation. These velocities should be used in the context 
of a full numerical model for planetesimal accretion that takes into account 
also the frequency with which different or equal size planetesimals collide. If 
the size distribution of the swarm privileges impacts between equal size plan- 
etesimals, it comes out clearly from the data in the tables that accretion is 
possible also for very close and eccentric binary systems. On the other hand, 
if collisions between different size planetesimals dominate the evolution of the 
planetesimal swarm, then the different alignment of perihelia might halt the 
accretion process. Unfortunately, the size distribution of a planetesimal swarm 
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is poorly constrained at the beginning of the accretion process and it may also 
be a misleading concept (Wetherill & Inaba, 2000). The frequently adopted 
assumption that "initial" planetesimals have a given size Rpi is numerically 
convenient but certainly an oversimplification. Regardless of the details of the 
planetesimal formation process, there have to be a dispersion ARpi in plan- 
etesimal sizes. The amplitude of ARpi is difficult to estimate, but it seems 
likely that it spans over at least one (if not several) order (s) of magnitude 
(see for instance Fig. 3 of Weidenschilling. 2000). If this is the case, then in 
the binary stars cases explored here, encounters leading to erosion ("NA" type 
encounters) might largely dominate over those leading to accretion (type "A"). 
We have checked this by performing a simplified test where we assume for the 
planetesimal size distribution a Gaussian function centered on ii = 5 km of 
variance ARp[. By counting the number of A-type encounters vs. NA-types 
using the values of (Av) (^jj^ ji^^ reported in Tabs.l and 2 we find that accretion 
dominates over erosion only when ARpi < 0.8 km. This value of ARpi appears 
to be small, possibly too small for any realistic initial size distribution. As a 
consequence, gas drag tends to prevent, or at least slow down, accretion in the 
initial stages of planetesimal evolution. As accretion proceeds, larger bodies 
are formed and impacts occur more and more frequently between diff'erent size 
bodies. The different phasing caused by gas drag for different size planetesi- 
mals increases the number of NA-type encounters and erosion dominates, at 
least for the binary parameters given in the tables. 

The limit of a purely dynamical approach that estimates the planetesimal 
relative velocity but neglects their size distribution can be partly overcome by 
taking into account that on average the projectiles delivering the maximum 
kinetic energy to a target of size i?2 are roughly those of size Ri ~ l/2i?2 (see 
Thebault et al., 2006). With this assumption we can probe the feasibility of 
accretion by testing (Av) against Vero only for a limited number of different 
size body pairs. By exploiting this simplification, Thebault et al. (2006) have 
tested the chances of planetesimal accretion for 120 different binary systems 
with semimajor axis Ub ranging from 10 to 50 AU and eccentricity from 0.05 to 
0.9. The mass ratio has been kept fixed to mt, = 0.5. In all the simulations only 
two size pairs have been considered, i?i = 2.5km vs. R2 = 5km, representative 
of small planetesimals, and i?i = 15km vs. R2 = 50km, representative of 
intermediate size bodies in the swarm. In Figs.5and6 different color coding 
are used to outline the planetesimal coUisional fate, depending on the binary 
parameters: 

• 1) Dark green. The evolution of a planetesimal disk in binaries belonging 
to these regions is not significantly affected by the companion star pertur- 
bations. {Av) remains low and accretion should proceed as in the standard 
model for single stars. 

• 2) Light green. Within these regions, {Av) is increased by the gas drag 
coupling to the secular perturbations but remains below Vgro- Accretion 
is still possible but the amplitude of the gravitational focusing factor 
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(^esc(fl2)/^'^)^ is significantly reduced. This slows down the accretion rate 
compared to the single star unperturbed case. 

• 3) Yellow. Here {Av) is comparable to the erosion velocities computed with 
the different prescriptions adopted at the beginning of the study for Vgro- 
This is a limiting situation where the preference of the system towards 
accretion or erosion depends on the details of the still not well defined 
collisional physics. We prefer not to draw any definite conclusion in this 
case. 

• 4) Orange and Red. In these regions, (Av) exceeds Vgro for all three dif- 
ferent collision outcome prescriptions considered. Collisions always result 
in a net mass loss and mutual accretion is impossible. 

As expected, systems with higher eccentricity or lower binary separation 
are more critical for planetesimal accretion. An empirical fit has been per- 
formed in Thebault et al. (2006) to outline the region in the (ab,ef,) phase 
space where accretion is possible for the two different size ranges of planctes- 
imals: 

66 ~ 0.013 (y^u) ' small 1 < R < 10km planctesimals, (10) 

and 

— - — ) , for large 10 < R < 50km planctesimals. (11) 

By extrapolating the fit to large values of one can figure out that for binary 
separations Ub > 90AU the planetesimal accretion process is not significantly 
perturbed by the companion star gravity. 

One interesting feature is the extent of the "light-green" region where 
planetesimal accretion is not halted but possibly slowed down since the grav- 
itational focusing factor is reduced. As a consequence, "classical" runaway 
growth is not possible and might be replaced by the alternative, possibly 
slower type II runaway mode identified by Kortenkamp et al. (2001). 



6 Type II runaway growth 

In this section we describe details of planetesimal accretion simulations in the 
Sun-Jupiter "binary" system (Kortenkamp & Wetherill, 2000a; Kortenkamp 
et al., 2001). This system has traditionally been used for initial characteriza- 
tions of dynamical effects in binary-star systems (e.g. Heppenheimer, 1978; 
Whitmire et al., 1998). In the models described here we use units of binary 
separation D rather than AU. Note also that the modeling included Saturn 
as well (a triple system) but we omit Saturn from most of the discussion, 
referring simply instead to the Sun-Jupiter "binary." 

Initial orbits for the planctesimals were circular, coplanar with the nebular 
midplane, and uniformly distributed from 10-50% of D. As the system evolves 
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Fig. 5. Encounter velocities averaged, over the time interval < t < 2 x 10'*yrs, 
between Ri = 2.5 and R2 = 5 km bodies at 1 AU from the primary star, for differ- 
ent values of the companion star semi-major axis and eccentricity. The short black 
vertical segments mark the limit beyond which < Av(iii,R2) > values correspond to 
eroding impacts for all tested collision outcome prescriptions. 

the combination of gas drag and secular perturbations from the companion 
Jupiter (and Saturn) leads to a size-dependence in e and / as well as in the 
phasing of the orbital orientation angles, as described earlier in this review 
chapter. Note that the system is fully three dimensional. The companion orbits 
are inclined with respect to each other and the protoplanetary disk. 

As described above, the distribution of encounter velocities is a critical fac- 
tor which determines how a population of planetesimals will accumulate into 
larger bodies. High encounter velocities can actually lead to impact erosion of 
the target rather than growth. On the other hand, low encounter velocities can 
gravitationally enhance the collision cross-section of the target (Opik, 1951). 
The effective cross-section becomes 

(Teff = a f 1 + , (12) 

where a is the target's geometric cross-section, Wcnc is the projectile's en- 
counter velocity, and Wosc is the combined surface escape velocity of the target 
and projectile. 
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Fig. 6. Same as Fig. 5, but for Ri = 15 and R2 = 50 km bodies. The short black 
vertical segments mark the limit beyond which < ^U{j;i,fl2) > values correspond 
to eroding impacts for all tested collision outcome prescriptions. The short dashed 
vertical segments mark the position beyond which orbital crossing occurs, despite 
the effects of gas drag, for 50 km bodies. 

The size-dependent phasing of orbital elements leads to low encounter 
velocities between similar size bodies and high encounter velocities between 
bodies of markedly different size. For the Sun-Jupiter model described above 
"binary" encounter velocities were calculated for all intersecting orbits over 
the entire range of planetesimal sizes. Figure 7 shows examples of these en- 
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Fig. 7. Mean encounter velocities (±cr) are shown 
for various size planetesimals with intersecting or- 
bits near 33% D, where D is the Sun- Jupiter sepa- 
ration distance. The size-dependent orbital phas- 
ing of bodies on intersecting orbits leads to high 
encounter velocities between bodies of markedly 
different sizes. The surface escape velocity of the 
combined target/projectile is shown by the hori- 
zontal dashed line. In the top two panels the ero- 
sion limit is indicated by the horisontal solid line. 
For encounter velocities above the erosion limit 
impacting projectiles crush and eject more than 
their own mass (crushing strength was assumed 
to be 10* ergs/g). Top: Encounter velocities for 
10^* g targets and projectiles of equal or lessor 
mass. Middle: Targets with masses of lO'^^ g are 
large enough that impacts of all projectiles lead to 
growth. Bottom: Planetsimals capable of growing 
to 10^** ! g have surface escape velocities near 500 
ms~^. Encounter velocities between these large 
planetesimals and bodies as much as 8-10 orders 
of magnitude less massive are below the escape 
velocity, allowing for significant gravitational fo- 
cusing. 

counter velocities for various target bodies and equal or lesser size projectile 
bodies with orbits near 33% D, where, again, D is the Sun- Jupiter separation 
distance. 

Modeling the growth of planetesimals in this scenario is considerably more 
difficult than tracking their orbital evolution. Full-scale iV-body simulations 
of planetesimal growth that include mutual gravitational perturbations, sec- 
ular perturbations, and gas drag arc beyond the reach of current techniques. 
Theoretically one would need to include ~10^^ small (^10^^ g) planetesimals 
to form a single 10^^ g terrestrial planet embryo. Direct A''-body integrations 
of mutually perturbing planetesimals cannot even remotely approach this fig- 
ure, treating only about ~10^ bodies over the time scale required. On the 
other hand, existing statistical simulations of planetary growth are not lim- 
ited by the number of bodies. However, these simulations assume that the 
orbits are completely randomized so they cannot easily accommodate the 
size-dependent orbital evolution described above. To overcome these short- 
comings, Kortenkamp et al. (2001) developed a hybrid approach that capi- 
talizes on the strengths of each technique. They used A^-body integration to 
map the size-dependent velocity distributions. These velocity distributions are 
then used in modified statistical accretion simulations to model the coUisional 
accumulation of the planetesimals. This approach led to the identification of 
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Fig. 8. Growth of planetesimals in a circumsolar disk perturbed by the "companion" 

Jupiter at a distance D from the sun. The habitable zone, simply defined here as 
the Venus-Mars region, stretches from about 10-25% of D. Initially all planetesimals 
have identical masses of 10^* grams. The initial surface density of planetesimals 
at 15% of D is 10 g cm~^, scaling as r~^^'^ with heliocentric distance r. This is 
roughly consistent with the so-called "minimum mass" surface density. Collision 
fragments smaller than 10 g (~1 meter) arc presumed lost via nebular gas drag 
and therefore removed from the simulation. Growth is calculated in four different 
regions centered on 15, 25, 33, and 40% of D. For comparison, the masses of Earth 
and Mars are iudic-ited by their respective symbols. The top two panels represent 
growth in the habitable zone and show initial orderly growth transitioning to Type 
II runaway growth of Mars-size planetary embryos in 10^ to 10® years. Beyond the 
habitable zone (bottom panels) growth is slower but eventually also produces Maxs- 
size embryos. Figure is adapted from Kortenkamp et al. (2001). 





Four separate simulations were performed, centered on 15%, 25%, 33%; 
and 40% of D. For comparison, the habitable zone extends from about 10- 
25% of D (crudely defined as the Venus-Mars region). In all four regions 
the growth is characterized as "orderly" (non-runaway, see Safronov, 1969) 
until the distribution reaches approximately lO-^** g, or about the size of the 
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largest asteroid — (1) Ceres. The planetesimal size distribution then becomes 
bi-modal, transitioning to runaway growth and producing Mars-size embryos. 
All planetesimals were assumed to be uniformly distributed across a region 
and any bodies separated by more than 10 mutual Hill radii were considered 
dynamically isolated (Chambers ct al., 1996) and not allowed to collide. This 
resulted in multiple runaway embryos emerging in each region. 

Note that mutual perturbations between planetesimals were not included 
in these calculations. This new form of runaway growth arises when secular 
perturbations and gas drag act together to establish size-dependent encounter 
velocities that remain low when colliding bodies are of similar size. Collisions 
between bodies of significantly different size arc at high velocity and can 
lead to cratering and erosion, but these simulations show that growth over- 
comes erosion (Fig. 8). This alternative mode of growth, which was eloquently 
dubbed "Type 11" runaway growth, allowed for the formation of terrestrial 
planet embryos throughout the habitable zone of the Sun- Jupiter "binary." 
This general result should apply regardless of whether the perturbations are 
from Jupiter-like companions formed earlier by disk-instability (e.g. Boss, 
1997; Mayer et al., 2002), stellar-mass objects in multiple-star systems, or 
perhaps even short-lived instabilities that lead to asymmetries in protoplane- 
tary disks. 

Classical "Type I" runaway growth occurs in a self-gravitating population 
of planetesimals. Random orbital kinetic energy is exchanged during gravita- 
tional encounters between large and small bodies and the population trends 
towards energy equipartition, a process dubbed "dynamical friction" (Stew- 
art & Kaula, 1980). Dynamical friction lowers the encounter velocities of the 
larger bodies with respect to each other, enhancing their effective collision 
cross-sections and increasing the rate at which they accumulate each other. 
Under these conditions nearly the entire growth period up to embryo-size is 
in Type I runaway mode. In our simulations, which arc not self-gravitating, 
the size-dependent phasing of orbital elements holds encounter velocities low 
between all similar-size bodies (typically 1 to 10 ms~^, see Fig. 7). Initially 
these encounter velocities exceed the planetesimal escape velocities so there 
is no enhancement of collision cross-sections and growth is orderly. As larger 
and larger bodies grow their escape velocities approach and then exceed the 
relatively low encounter velocities, causing the transition from orderly growth 
to Type II runaway growth. In this way, the effects of dynamical friction are 
mimicked by the size-dependent phasing of orbital elements. 

7 Conclusions 

The process of planetary formation around the primary star of a binary sys- 
tem is complicated in all its stages by the gravitational perturbations of the 
companion star. However, the c;xistence of some gas giant planets in binary 
systems with separation of a few tens of AU suggests that these perturba- 
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tions are not always strong enough to prevent the formation of a planet. We 
have analysed within the standard model of planet growth how the seciilar 
perturbations of the companion star affect the various stages of planetesi- 
mal accretion. Ours studies show that when we include also the effects of gas 
drag on the motion of small planctesimals, their orbits become phased both 
in eccentricity and perihelion longitude (Marzari & Scholl, 2000). At a first 
sight, this alignment favors fast accretion keeping low the encounter velocities 
between the bodies in spite of the large values of eccentricity forced by the 
secular perturbations. However, the angle towards which the perihelia align 
depends on the size of the planctesimals and on their distance from the star 
(Thebault et al., 2004, 2006). In the earlier stages of planctesimal accretion all 
the bodies have roughly the same size and the encounter velocities, because 
of the favorable phasing, are low. However, as soon as larger bodies emerge in 
the population the different alignments of the orbits lead to large impact ve- 
locities that may prevent the onset of runaway growth or even cause erosion of 
the bodies inhibiting the formation of a planet. Through a full numerical ap- 
proach we have analysed to which extent the different alignment influences the 
random planctesimal velocities. We have also mapped the values of the binary 
orbital parameters (semi-major axis and eccenticity) for which accumulation 
dominates over erosion despite of the secular perturbations, or viceversa. Our 
modelling shows that for binary of separations Ub < 40AU, only very low Cb 
allow planctesimal accretion to proceed as in the standard single-star case. 
On the contrary, only relatively high values (of at least 0.2 in the closest 
ttb = lOAU sepration explored and at least ~ 0.7 for Ub = 40 AU) lead to 
a complete stop of planctesimal accretion. In most cases, when the pertur- 
bations of the massive companion on the planctesimal disk is significant but 
not strong enough to halt accretion, runaway growth can still occur, but in a 
different way respect to the classical "Type I" runaway typical of planctesimal 
populations around single stars. This new type of growth termed "Type H" 
runaway allows planet formation to occur in binary-star systems with much 
tighter orbits than previously suggested (Heppenheimer, 1978; Whitmire et 
al., 1998). However, there is (at least) one caveat. As already noted, the growth 
simulations described above and represented in Figure 8 did not include mu- 
tual perturbations of the planctesimals themselves. Crude attempts at includ- 
ing self-gravitating planctesimals (Kortenkamp & Wetherill, 2000b) indicated 
that when the size distribution reaches lO'^'* to 10^^ g the mutual perturba- 
tions are beginning to become important, although they are still dominated by 
the size-dependent phasing of secular perturbations. This suggests that just 
as Type II runaway growth is getting underway accretion may either stall or 
perhaps transition to the classical Type I runaway growth. To explore these 
possibilities we are working to modify the multi-zone planctesimal accretion 
code of Weidenschilling et al. (1997) to include secular perturbations from 
massive companions as well as mutual perturbations from planctesimals and 
nebular gas drag (Kortenkamp et al., 2006). 
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